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Abstract. As additional perturbative degrees of freedom, it is known that magnetic fields of 
inflationary origin can source curvature perturbations on super-Hubble scales. By requiring 
the magnetic generated curvature to remain smaller than its inflationary adiabatic counter- 
part during inflation and reheating, we derive new constraints on the maximal field value 
today, the reheating energy scale and its equation of state parameter. These bounds end up 
being stronger by a few order of magnitude than those associated with a possible backreaction 
of the magnetic field onto the background. Our results are readily applicable to any slow-roll 
single field inflationary models and any magnetic field having its energy density scaling as a 1 
during inflation. As an illustrative example, massive inflation is found to remain compatible 
with a magnetic field today Bq = 5 x 10~ 15 G for some values of 7 only if a matter dominated 
reheating takes place at energies larger than 10 5 GeV. Conversely, assuming 7 = —1, massive 
inflation followed by a matter dominated reheating cannot explain large scale magnetic fields 
larger than 10 -20 G today. 
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1 Introduction 

Current measurements of magnetic fields filling the intergalactic medium show that they are 
not vanishing on the largest length scales. Using the spectra of blazars, and under the safest 
assumptions, Ref. [1] finds the two-sigma limit Bq > 10 -17 G whereas other data yield that 
Bq > 5 x 10~ 15 G [2-4] and Bq > 10~ 20 G [5]. A non- vanishing magnetic field today, and 
on the largest length scales, suggests that it has a primordial origin. Indeed, as any other 
perturbation modes, large scale magnetic fields are necessarily of super-Hubble wavelengths at 
higher redshifts, and this triggers the problem of an astrophysical generation mechanism [6-8]. 
Since inflation solves exactly the same problem for the generation of the primordial curvature 
perturbations, it sounds the best candidate to generate magnetic field [9]. However, in flat 
space [10], the conformal invariance of electromagnetism prevents an amplification mechanism 
to take place during inflation. Even worse, as the magnetic energy density redshifts as 
p B oc 1/a 4 , its contribution during inflation can become dominant and leads to a severe 
backreaction problem [11-15]. In fact, backreaction can already appear during the reheating 
era as soon as the mean equation of state parameter u; r ch < 1/3. As shown in Ref. [16], 
this yields some non-trivial constraints between the energy density p re h at the end of the 
reheating era, the equation of state parameter w r< ^ and the present Hubble scale value of 
the magnetic field Bq. As a result, if magnetic fields have an inflationary origin, conformal 
invariance has to be broken, at least during inflation [17-27]. 

Assuming the background evolution is under control, one still has to consider the gravi- 
tational effects of the magnetic degrees of freedom onto the evolution of the cosmological per- 
turbations. Among others, one expects the generation of primordial non-Gaussianities [28-32] 
and a non-conservation of the curvature perturbation on super-Hubble scales, in a way simi- 
lar to the presence of entropy modes [28, 33-35]. In this context, magnetic fields should not 
induce too large curvature perturbations on super-Hubble scales which would otherwise spoil 
the standard adiabatic contribution of inflationary origin. In this paper, we look into this 
issue, both during inflation and reheating, and use a phenomenological model to describe the 
evolution of the magnetic field during inflation as B oc a 7 . Although our approach does not 
allow to derive the backreaction of the perturbations onto the magnetic field itself, our results 
should give the correct order of magnitude. In particular, and in addition to 7, we show that 
the maximal allowed value of Bq on Hubble scales today depends on the way the reheating 
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proceeded and on the inflation model itself. By assuming only slow-roll during inflation and 
a constant equation of state during reheating, we derive a generic formulae giving the upper 
bound of 5o hi terms of the energy scale of reheating p l J e ^ and the corresponding equation of 
state parameter w re h- Such a result therefore generalizes the background constraints derived 
in Ref. [16] and extend the results of Ref. [34] to any reheating history. Our findings are then 
applied to large field massive inflation and the small field inflationary models. For instance, 
we find that massive inflation can only be compatible with a matter dominated reheating 
and a magnetic field value today of Bq = 5 x 10~ 15 G if the energy scale of reheating is higher 
than -EVch ^ 10 5 GeV, independently of 7. If some assumptions are made on 7, the bounds 
can be much stronger and are represented in Fig. 1. Our main formulae are Eq. (2.7) and 
Eqs. (2.11) to (2.13) which explicitly give the amplitude of the curvature perturbation given 
an inflationary potential, the magnetic field today and the reheating parameter. 



2 Super-Hubble curvature perturbations from magnetic fields 

The super-Hubble scale curvature perturbation £ on the constant energy density hypersurface 
sourced by an inhomogeneous primordial magnetic field is given by [34], 

where we have imposed an initial condition that C(t*) = 0- I n the actual situation, i* may 
be taken to be a Hubble crossing time. The first term, which we denote by £1, represents the 
contribution due to the relative entropy perturbation, 5P rc \ = 6Pb — ^5pb where 5Pb and dps 
denote pressure and energy density perturbation of the electromagnetic field. The quantities 
P and p are the total pressure and the total energy density, respectively. Assuming that 
the electromagnetic field obeys the Maxwell equations after the end of inflation, £1 becomes 
independent of time after reheating. On the other hand, the second term, which we denote 
by £2 , is sourced by the anisotropic stress of the magnetic field. Since the anisotropic stress 
persists even after reheating, £2 still continues to evolve during and after reheating. 

2.1 Phenomenological model 

In order to evaluate Eq. (2.1), we need to specify the time evolution of 5P re \ and II during 
inflation, which requires specification of the generation model of the magnetic field. Since we 
do not want to concentrate on the particular model of the magnetogenesis, in this paper, we 
will take a phenomenological approach and make the following ansatz: 



5P rel {t) = ap Bcnd 



a(*) 

^cnd 



7 



n(t) = Ppe 



^cnd 



(2.2) 



where a, (3 and 7 are constant parameters and pB end is the energy density of the magnetic 
field at the end of inflation, which is treated as first order perturbation. 

Then, adopting the slow-roll approximation and assuming a constant equation of state 
during reheating, Eq. (2.1) can be explicitly integrated in terms of the number of e- folds TV: 

/•-Wend p l{N-N end ) 1 „ rN teh /-i _q \ 

Ci = -ap B , / : dN--^ [i 6Wveh) e (-i+3- rch )(iV-^n d ) dNj 

JN, 4> 2 3 /)end JN Bnd 1 + ^reh 

(2.3) 
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where a dot stands for differentiation with respect to N. The first term can be further 
simplified by assuming slow-roll while the second can be integrated explicitly. In reduced 
Planck mass units (M^ = 8nG = 1), and at first order in the Hubble flow functions, one gets 



Ci = 3ap E 



In this expression 
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(2.4) 



(2.5) 



^reh \ Prch / \ Pend , 

is the reheating parameter [36-39]. The labels "end" and "reh" denote respectively the end 
of inflation and the end of reheating, i.e. the beginning of the radiation era. Since after 
inflation we assume the magnetic field to obey Maxwell equations, one can express the ratio 
pB ead /Pend m terms of energy densities today [16]. In particular, one has Ps cnd = Pb (1+^end) 
where z en d is the redshift at the end of inflation. Assuming instantaneous transitions [40] 
and making use of the reheating parameter, one has 
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(2.6) 



<?0 /3 3reh/ (^/h 9o) IS the 



Here p yo = 3Hq£1® is the total radiation density today, and Q r ch 
measure of the change of relativistic degrees of freedom between the reheating epoch and 
today, where q and g respectively denotes the number of entropy and energetic relativistic 
degrees of freedom at the epoch of interest. Plugging Eq. (2.6) into Eq. (2.4), and making 
use of the Friedmann-Lemaitre equations to express p nd in terms of the field potential 
V cn d = V(4> cnd ), one finally gets 
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(2.7) 



Here AN((f>) < is the number of e-folds before the end of inflation and can be obtained 
from the slow-roll trajectory 



A7V(0) = N(<f>) - N ( 



end 



ad V ,f 



(2.« 



The first Hubble flow function ei(<p) is evaluated along the field trajectory and is also uniquely 
determined by the potential in the slow-roll approximation [41] 



1 fv. 



2H 2 2 V V 



(2.9) 



By definition, the quantity ei cn d = ei(<^end) is unity for inflationary models ending by slow- 
roll violation but can be much smaller than unity for inflationary models ending by tachyonic 
instabilities. This formula allows us to evaluate (j once the inflation model and the reheating 
are specified. 
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In a similar manner, straightforward calculations yield £2 at any e-fold N during the 
radiation era. For this, it is convenient to split the integral over the three domains, inflation 
for N* < N < N cn a, reheating for iV en( j < N < N rc ^ and radiation era with N > N^: 



UN) 



c (inf) +c (reh) +c (rad) w 



(2.10) 



Again assuming slow-roll during inflation, and a constant equation of state during reheating, 
one gets, at first order in the Hubble flow functions: 
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(2.13) 



As before, €i(4>) stands for the first Hubble flow function evaluated along the field trajectory 
and ei cn d = ei(^ e nd)- The quantity AiV* = N* — N en ^ < is the number of e-folds before 
the end of inflation at which the pivot scale k* crossed the Hubble radius during inflation. It 
does not depend on in the previous equations but, as discussed below, it is an inflationary 
model-dependent function of the reheating parameter -R m d (see Ref. [39]). The same remark 
holds for V* = V(<fi*). Finally, when relevant, we have dropped all terms involving ei* as we 
always have ei* <C 1 for all inflationary models. 

Let us stress that the quantity Q r (N) grows during the radiation era, as N — A^ en . 
However, Eq. (2.13) only makes sense if the perturbation mode under consideration is super- 
Hubble. As a result, an upper bound for N — 7V ro h is given by the number of e-fold after the 
end of reheating at which the pivot scale k* re-enters the Hubble radius. 

Finally, let us notice that the expressions for £1, C2 an d are an propor- 

tional to the factor p B ,/p e nd — R>7adP B al P~ia ■ As shown in Ref. [16], imposing this factor to 
be smaller than unity is equivalent to avoid magnetic field backreaction over the background 
energy density during reheating. As here we are requiring that |£| < 10~ 5 , our constraints 
are expected to be typically, and at least, 2.5 order of magnitude stronger than those derived 
from the background evolution (see Fig. 1). 



2.2 Reheating consistent slow-roll 

As already mentioned, both (1 and £2 depends on parameters of the inflationary model as 
well as the reheating energy scale. In particular, one needs to determine the value of all "*" 
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quantities, i.e. evaluated AiV* e-folds before the end of inflation. These parameters cannot 
be freely chosen as one wants them to be compatible with the observed amplitude of the 
CMB fluctuations. In fact, as shown in Ref. [36], AN* is itself a function of -R ra d- I n order 
to ensure consistency with reheating, we here follow the slow-roll approach of Ref. [39] that 
we briefly recap. 

The e-fold iV* is by definition solution of k*/a(N*) = H(N*) during inflation. This 
equation can be recast in terms of "observable" quantities, namely 



— (l + ^end)e 

a 



(2.14) 



The right hand side of this equation can be fixed by the amplitude of the adiabatic primordial 
power spectrum, which is a well measured quantity. At leading order in slow-roll 
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From Eqs. (2.5) and (2.6), 

Zend can be expressed in terms of i? ra d an d Pend; the energy 
density at the end of inflation. The latter can, in turn, be further simplified using the 
Friedmann-Lemaitre equations for a scalar field (in Planck units) 
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(2.16) 



As before, one can drop the term in ei* <C 3. The advantage of this last expression is that 
it does no longer depend on the potential normalization but involves only H*. Plugging 
everything back into Eq. (2.14), AiV* is a solution of the algebraic equation [36] 



AiV* = -lni? rad + iVo + j hi 



Knd 



ei*(3 — ei end ) V* 



1 



ln(87r 2 P, 



(2.17) 



Here the quantity |iVo| roughly measures the number of e-folds of deceleration and is defined 
in Planck units by 

k*/a 
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(2.18) 



Notice that the trajectory is needed to evaluate the right hand side of Eq. (2.17) as V* and 
ei* are functions of 0(iV*). One can nevertheless render this equation more explicit in terms 
of (/>* and prch by assuming slow-roll and expanding i? ra( j from Eqs. (2.5) and (2.16), i.e. 
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Once the slow-roll trajectory is integrated, this expression is explicit in cp* and can be solved 
by specifying only E rc h = p l JJ^ and w TC ^. Plugging the result into Eqs. (2.11), (2.12), (2.13), 
(2.19), (2.20) and finally Eq. (2.7) allows to determine Q{N*)/ p BQ uniquely from the input of 
E^ch and w re h- Imposing that |C(iV*)| 2 <g; actually yields the reheating-dependent upper 
bound on p B . 



3 Application to some representative models 
3.1 Large field models 

The potential energy for the large field models is given by 

V(<f>) = MV, (3.1) 

where M is a constant of mass dimension (in Planck unit) and p is a positive number. In order 
to be definite, we assume that this potential is correct not only for large <f> responsible for 
inflation but also for small (f> relevant for oscillating period and reheating. For this potential, 
the field value at which inflation terminates is given by 

<?Wd = , (3.2) 

the solution of ei(</> en d) = 1 where 

eiWO-l^- (3-3) 

The slow-roll evolution of <f> during inflation is given by integrating Eq. (2.8). In Planck 
units, one gets 

AN = kp ( ^ nd " ^ ' (3 ' 4) 

where, as before AN = N — N en ^, is the number of e-fold measured from the end of inflation. 
After inflation, <fi oscillates around the minimum such that the natural equation of state 
parameter is given by w TC h = (jp — 2)/(p + 2) [42]. 

Plugging Eqs. (3.1), (3.2), (3.3) and (3.4) into Eq. (2.7) and Eqs. (2.11) to (2.13) gives 
the expression of £(A*)/p B() in terms of the potential parameter p, i? ra d and 0* only (plus 
the magnetic parameters). What remains to do is to solve Eq. (2.20) to get in terms of p 
and /9 ro h, from which -R r ad is determined by Eq. (2.19). 

For the large field model, an analytic solution of Eq. (2.20) can be found in terms of 
the Lambert function [39], but as the various integrations entering into the expression of 
((N*) that can only be performed numerically, we have here preferred to solve this equation 
numerically. The various cosmological parameters have been set to their preferred values 
from the WMAP data [43-45], i.e. P* ~ 2.16 x 10~ 9 , h ~ 0.72, 0° = 4.6 x 10~ 5 , and we 
have set all Q to unity for simplicity (they have only a small effect). With a pivot scale 
chosen at k* = 0.05 Mpc" 1 , one has iVo — —62. The final result is a parametric curve 

#0 = #a,/3, 7 ,p(-Erch,W roh ), Solution of 

C(N*) = VK, (3.5) 

which separates the plane (-E? re h, -Bo) in two regions. Above this curve, the combination of 
the magnetic field value today and the energy scale of reheating would be such that the 
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Figure 1. Limits on the amplitude of the present magnetic field of inflationary origin with respect 
to the reheating energy scale. We have assumed a large field model with p = 2 (and generic values for 
a = 4/3, /3 = 1). Each line is the solution of Eq. (3.5) for various values of 7. The region above the 
corresponding line is excluded. The reheating energy scale E rc h = p 1 ^ varies from 100 McV (BBN 

energy scale) to the energy at which inflation ends, see Eq. (2.16). For convenience, the CMB 
lower bound for the large field with p = 2 reheating energy has been reported £? re h > 70GeV (95%), 
see Ref. [38]. We have also represented the region excluded by magnetic field backreaction over the 
background energy density during reheating, see Ref. [16]. 



magnetically generated super-Hubble curvature perturbation would be equal or larger than 
the adiabatic modes generated during inflation. The allowed region therefore lies under this 
curve. Let us notice that this is a very conservative upper limit as current constraints on 
isocurvature modes show that they cannot exceed 10% of the adiabatic counterparts [46, 47]. 

In Fig. 1, we have plotted these limits for the large field model with p = 2 and for 
various values of 7 (taking a = 4/3 and f3 = 1 as reference values of 0(1) parameters). 
Negative values of 7 provides the tightest bounds on the magnetic field. This is expected 
since in that situation the magnetic contribution increases more and more for a — > 0, i.e. 
deep during inflation. As it is already known, standard values for 7 = —4 during inflation 
would even generate a strong backreaction problem [11, 12, 15]. Here, even for 7 = — 1, 
curvature perturbations generated by the magnetic fields become important during inflation 
and C(-W*) — Ci + (2 • We see on the figure that, for 7 = — 1, one cannot actually generate 
a magnetic field today larger than Bo = 10~ 20 G. In the opposite situation, 7 > 0, deep 
during inflation the magnetic effects are very small and ((N*) ~ + Ci™^ ( N *) is mostly 

generated after inflation. Remembering that after inflation p B oc a~ 4 , at fixed Bq, a longer 
period of reheating, i.e. a lower reheating temperature, results in higher magnetic field 
values at the end of inflation; and hence a larger £. As for the background case discussed in 
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Figure 2. Limits on the amplitude of the present magnetic field of inflationary origin with respect 
to the reheating energy scale for a large field model p = 2, 7 = 2 and for various values of w rc h- 



Ref. [16], this effect is enhanced if the energy density of the universe decreases more slowly 
than radiation, as this is the case for w re h = here. As can be checked in Eqs. (2.12) and 
(2.13), the C-de-pendency of these two terms with respect to 7 is very weak (as opposed to 
C^ nf) ). In Fig. 1, one indeed sees that the bounds become almost insensitive to 7 as soon as 
7 > 2 such that this limit is actually very conservative and somehow model-independent. At 
around -E re h = 10 14 GeV, a spiky bump is observed for 7 = 1 curve (and similar one for 7 = 
curve at higher value of E rc h.). This is due to the occasional vanishing of £ ; which is possible 
since £ is a function of i? ra d an d hence of E rc h. At that point, £ vanishes irrespective of the 
amplitude of the magnetic field and, as a result, we have no constraint on the magnetic field, 
which shows up as the bump. 

In Fig. 2, we have plotted the magneto-reheating constraints, still for a large field 
model with p = 2, but assuming another equation of state parameter tt> re h- As expected, the 
more w re h becomes negative, the tighter the limits are whereas there is no constraints for a 
radiation-like reheating era. 

3.2 Small field models 

Let us next consider the small field models for which the potential is given by 



V(4>) = M 4 



1 



(3.6) 



where M and /x are constants of mass dimension in Planck unit and p is a positive number. 
Notice that this potential is valid only for the inflationary regime (p/fM<^.l and has to be 
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p/M -=10 p=2 w re =0 




Figure 3. Limits on the amplitude of the present magnetic field of inflationary origin assuming the 
small field models with p = 2 and fi = 10M P! . For each labelled value of 7, the region above the 
corresponding line is excluded. These bounds are almost unchanged for larger values of \i (see text). 



replaced by another function after the end of inflation. If this is a quadratic function, one 
can assume u> rc h — 0. For this potential, normalizing the field value by \x as x = the 
slow-roll trajectory reads (in Planck units) 
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where Xend is determined by ei(xend) = 1- From 
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one finds a transcendental equation for x 
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Equation (3.7) can also be applied to the case p = 2 by taking the limit p — > 2, which yields 
terms in ln(x/Xend)- Along the lines detailed for the large field models, these expressions are 
enough to completely determine the upper bound Bq = B(E re h, w re h) parametrized by /j, and 
p (see Sect. 3.1). 

In Fig. 3, we have plotted the limits on the amplitude of the present magnetic field of 
inflationary origin as a function of the reheating temperature assuming small field models 
with p = 2 and \x = 10M Pi . Notice that for u) ren = 0, there is currently no bounds on the 
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Figure 4. This panel shows the limit on the amplitude of the present magnetic field of inflationary 
origin assuming the small field models with with various values of p and at fixed /i = M p! . As before, 
the region above the lines is excluded and a = |, /3 = 1. 



reheating energy scale coming from CMB alone [38]. Let us also recap that the models having 
p = 2 are compatible with the current spectral index and tensor-to-scalar ratio only for large 
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values of //, typically // 3> M p; [36]. The bounds of Fig. 3 exhibits the same behaviour as the 
large field models, i.e. tighter bounds on the magnetic field for lower reheating temperatures 
and for lower values of 7. We also find that those results remain insensitive to larger values of 
//. This is not surprising since, for the small field models, the first two-Hubble flow functions 
become independent of /i in the large //-limit [38]. Moreover, the above equations imply that 
1 — Xend = 0{M pl /fj,) and Xend — X* = ®(M pl / p). Therefore the integral in Eq. (2.7) scales 
as Mp //i 3 and hence £i becomes also independent of // in the large //-limit. 

In the top panel of Fig. 4, the limits on the amplitude of the magnetic field have 
been derived for the case p = 5, and we have made both 7 and // vary. Contrary to the 
case p = 2, sub-Planckian values for /t can be made compatible with CMB data [38]. The 
magneto-reheating bounds become generically stronger at small //; but such a behaviour can 
be transiently inverted for some very large values of 7 3> 1. For instance, this is the case for 
7 = 50, but for values of \i < W~ 2 the upper bound moves again downwards in Fig. 4 (not 
represented). As before, for 7 < the bounds are driven by the behaviour of £1 because the 
magnetic effects are dominant deep in inflation. For 7 > 0, magnetic effects during reheating 
are the most important, and Q r * s ^ ne dominant term such that the //-dependence ends 
up being related to the values of 7. In the limit fj,/M Pl <C 1, one can nevertheless use some 
crude approximations to guess the dependency in //. One has x* — \p(p~ 2)// 2 /|AiV*|] 1 /( p-2 ) 
such that the dominant terms in Eq. (2.7) scale as (in Planck units) 

dx oc //- 2 p/(p- 2 ) . (3.10) 

Similarly, one can use Eq. (2.19) to get i^ ad oc ei <^ch-i)/(3+3™ roh )^ Jn ^ Hmit ^ < ^ 

Eq. (3.8) shows that ex* oc // 2 p/(p~ 2 ) and, for w rc h = 0, one finally gets that £1 oc ex* 2 ^ 3 - As 
expected, this quantity increases as // decreases. For 7 > 0, the behaviour of £2 is now driven 
by powers of R ra d, which again increases when // decreases, albeit in a different way. 

In the bottom panel of Fig. 4, we have fixed // = M pl and plotted the upper limits for 
p = 3, p = 8 and p = 20. Since p only appears in the equation as an order unity factor, the 
final dependency in p remains weak. Notice again the change of behaviour between large and 
low values of 7 that for /t, be traced back to which part of £ contributes the most. 

4 Conclusion 

Recent observations of the cosmic rays reveal that the magnetic fields are ubiquitous in the 
universe. Lack of the convincing astrophysical explanation for the origin of such magnetic 
fields has led some theorists to seriously consider the possibility that those magnetic fields 
are produced during the primordial inflation. In addition to the obvious condition on the 
model of the inflationary magnetogenesis that it must produce the observed amplitude of 
the magnetic field at the relevant scales, there is another condition that must be taken 
into account for whatever the model of the magnetogenesis is. Since the inhomogeneous 
magnetic fields enter the right hand side of the Einstein equation as perturbations of the 
energy- momentum tensor, they induce the metric perturbation in addition to the standard 
adiabatic one produced by the inflaton fluctuation. The condition must be satisfied that 
this additional metric perturbation should not exceed the observed one in order for any 
inflationary magnetogenesis model to be embedded consistently in the standard model of the 
early universe. 



Xend i/4 

x* P 



p(p - 2) 



X 
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In this paper, we first derived general expression of the resulting super-horizon scale 
curvature perturbation sourced by the primordial magnetic field, without resorting to the 
specific model of inflationary magnetogenesis but with a phenomenological assumption that 
the magnetic field energy density and the associated anisotropic stress scale as a 7 , where 7 is 
a free parameter, so that our result provides a wide coverage for any magnetogenesis model 
satisfying the above scaling behavior. Our result can be also applied to any canonical single 
field inflation model followed by the oscillations of the inflaton with any equation of state and 
by reheating with any reheating energy scale. Given the inflation model and the measure- 
ment of the curvature power spectrum at the pivot scale, our formula allows us to put bound 
on the combination of the amplitude of the today's magnetic field, reheating energy scale 
and its equation of state parameter. In practice, our formula requires numerical solution of 
the algebraic equation to find the inflaton field value corresponding to the Hubble crossing of 
the pivot scale and one dimensional numerical integration for evaluating the curvature per- 
turbation, both of which are quite feasible to implement. Our perturbation bound is tighter 
by a few orders of magnitude than the one given in the literature, which is derived from 
the requirement that the magnetic field energy density be smaller than the background en- 
ergy density so that it does not destroy the homogeneous and isotropic Friedmann-Lemaitre 
metric. 

We then applied our formula to the large field inflation models and the small field 
inflation models, which are representative models of inflation. In either case, the upper bound 
on the value of the magnetic field today becomes tighter for negative value of 7, where the 
magnetic curvature perturbation is dominantly produced deep inside inflation. In particular, 
the magnetic field strength of 10~ 15 G, which has the observational relevance, with 7 < — 1 is 
completely incompatible with the perturbation bound. On the other hand, the upper bound 
is insensitive for positive 7 and saturates at 7 = 0(1). Thus we can interpret the upper 
bound for 7 = 0(1) as the fairly conservative bound. With 7 being fixed, the upper bound 
on the magnetic field becomes severer as the reheating energy scale becomes lower if the 
reheating equation of state parameter is less than 1/3. This is expected since the magnetic 
field energy density contributes more to the total energy density as we go back into the past. 
For example, for the chaotic inflation model with a quadratic potential, the magnetic field 
today must be less than 10 17 G if the reheating energy scale is as low as 70 GeV that is the 
possible minimal energy scale allowed by WMAP 7yr data. If the cosmic magnetic fields 
actually turn out to be of the inflationary origin in the future and search for the realistic 
model of the inflationary magnetogenesis becomes an indispensable pillar in constructing a 
realistic evolutionary scenario of the early universe, our perturbation bound derived in this 
paper provides one of the necessary conditions that must be considered for constraining the 
magnetogenesis models. 
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